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Abstract 

This article deals with invariant manifolds for infinite dimensional ran- 
dom dynamical systems with different time scales. Such a random system 
is generated by a coupled system of fast-slow stochastic evolutionary equa- 
tions. Under suitable conditions, it is proved that an exponentially tracking 
random invariant manifold exists, eliminating the fast motion for this cou- 
pled system. It is further shown that if the scaling parameter tends to zero, 
the invariant manifold tends to a slow manifold which captures long time 
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dynamics. As examples the results are applied to a few systems of cou- 
pled parabolic-hyperbolic partial differential equations, coupled parabolic 
partial differential-ordinary differential equations, and coupled hyperbolic- 
hyperbolic partial differential equations. 

Key Words: Stochastic partial differential equations (SPDEs); random 
dynamical systems; multiscale systems; random invariant manifolds; slow 
manifolds; exponential tracking property 
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1 Introduction 

The theory of invariant manifolds serves as a tool for analyzing dynamical behav- 
iors of deterministic systems. It was first introduced by Hadamard [20], then by 
Lyapunov [23] and Perron [29] for deterministic systems. It has been further de- 
veloped by many authors for infinite dimensional deterministic systems; see, e.g., 
[30] [21 [9l [121 [21] • More recently, invariant manifolds have been investigated for 
infinite dimensional stochastic systems; see [T71 [HI [231 [23 [3l [6l [151 [H] among 
others. 

Some systems evolve on fast and slow time scales, and may thus be modeled by 
coupled singularly perturbed stochastic ordinary or partial differential equations 
(SDEs or SPDEs). For SDEs with two time scales, SchmalfuB and Schneider [32] 
have recently investigated random inertial manifolds that eliminate the fast vari- 
ables, by a fixed point technique based on a random graph transformation. They 
show that the inertial manifold tends to another so-called slow manifold as the 
scaling parameter goes to zero. Qualitative analysis for the behavior of the slow 
manifold for slow-fast SDEs on the long time scales can be found in Wang and 
Roberts [33]. 

In the present paper, we consider invariant manifolds for stochastic fast-slow 
systems in infinite dimension. Namely we investigate the following system of fast- 
slow stochastic evolutionary equations, which could be coupled SPDEs, or coupled 
SPDEs-SDEs, 

= -Ax'' + - fix", y") + -^w, in Hi, 
e e y/e 

y' = By' + g{x',y'), in H2, 
where A and B are generators of Cq— semigroups, the interaction functions / and 

m 

g are continuous. The noise process 10=^1 ^j'^j^ where {w^jjjLi are two-sided 
Wiener processes (or Brownian motions) taking values in M and hj {1 < j < m) 
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are given elements in Hi. The small parameter e > representing the ratio of the 
two time scales. The precise conditions on these quantities will be given in Section 
131 which allow our framework to deal with multiscale coupled parabolic-hyperbolic 
systems and coupled hyperbolic-hyperbolic systems. 

It is worthy mentioning that in the situation we consider here, the noise pertur- 
bation of the fast motion equation is additive type. The reason is that the problem 
for existence of random dynamical systems generated by stochastic partial differen- 
tial equations with general multiplicative noise is still unsolved (for details see [T7]). 
The main goal in this paper is to establish, for e > small enough, the existence 
of a random invariant manifold with an exponential tracking property for the 
above stochastic system. Thus as a consequence, this system can be reduced to an 
evolutionary equation with a modified nonlinear term, which is useful for describ- 
ing the long time behavior of the original coupled stochastic system. There are 
usually two approaches to construction of invariant manifolds: Hadamard graph 
transform method (see [31], [T7]) and Lyapunov-Perron method (see [HI [181 [6] ) . We 
achieve our results by the latter which is different from the method of random graph 
transformation in [52]. In this approach one key assumption is that the Lipschitz 
constant of the nonlinear term in fast component is small enough comparing with 
the decay rate of the linear operator A. In particular, under suitable conditions 
it is further shown that this manifold can be asymptotically approximated for 
e sufficiently small by a slow manifold M° for a reduced stochastic system. We 
note that, in the case of Lyapunov-Perron method applied to a coupled stochastic 
systems, the existence of an random invariant manifold for the coupled stochastic 
parabolic-hyperbolic equations, that do not contain two widely separated char- 
acteristic timescales, is obtained by Caraballo Chueshov and Langa in [6]. We 
remark that, whereas the existence of slow manifold is not studied, in their paper 
the just mentioned authors also verify that this random manifold converges to its 
deterministic counterpart when the intensity of noise tends to zero. 

This paper is organized as follows. In Section 2, some basic concepts in random 
dynamical systems and random invariant manifolds are recalled. Our framework 
is presented in Section 3. In Section 4, we establish the existence of a random 
invariant manifold possessing an exponential tracking property, and then in 
Section 5 we show converges to a slow manifold M° with rate of order 1. 
Section 6 is devoted to a few illustrative examples. Remarks on local manifolds for 
systems with local Lipschitz nonlinearities are given in Section 7. 

2 Preliminaries on random dynamical systems 

We now recall basic concepts in random dynamical systems (RDS) and random 
invariant manifolds (RIM). For more details, see [Tl \T7\ [18]. 
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Definition 2.1. Let ^,P) be a probability space and 9 = {6t}t£R be a flow on 
Q which is defined as a mapping 

and satisfies 

• 9o = idn, 

• 9s9t = 6s+t for all s, t G M, 

• the mapping {t, uj) H- dtuj is (=^^(R) x — measurable and ^^P = P for all 
t G M. Then the quadruple {Q, J-", P, 6) is called a driving dynamical system. 

We will work on the driving dynamical system represented by Wiener process. 
To be more precise, let Q = Cq(R,W^) be the continuous paths u{t) on M with 
values M"* such that uj{0) =0. This set is equipped with the compact-open topol- 
ogy. Let J-" be the associated Borel a— field and P be the Wiener measure. Then 
we identify u with 

{wiit),W2it), ■ ■ ■ , Wmit)) = Uit), t G M. 

The operators 6t forming the fiow are given by the Wiener shift: 

etuj{-) = u{- + t)- u{t), u en, t G M. 

Note that the measure P is invariant with respect to the above fiow and then the 
quadruple {Q, J-", P, 6) is a driving dynamical system. 

Definition 2.2. Let (H, de) be a metric space with Borel a— field ^(H). A cocycle 
is a mapping: 

: M+ X X H ^ e, 
which is (^(]R+) x ^ x ^{M),^(B.)) -measurable such that 

0(0, u, x) = X, 

(pit + S, CJ, x) = (j)(t, 9sUJ, (f){s, u, x)), 

for t, s G M.^, u E Q, and a; G H. Then (p together with the driving system 9 forms 
a random dynamical system (RDS). 

A RDS is called continuous (difierentiable) if a; — )■ (j){t,u,x) is continuous (dif- 
ferentiable) for t > and w G fi. A family of nonempty closed sets M = {M{oj)} 
contained in a metric space (H, || ■ ||e) is called a random set if for every ?/ G EI the 
mapping 

w inf ||a;-i/||e 
is a random variable. Now we introduce the random invariant manifold concept. 
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Definition 2.3. A random set M{uj) is called a positively invariant set if 

(pit, u, M{u)) C M{etuj), fort>0, u e n. 
If M can be represented as a graph of a Lipschitz mapping 

ij{-,u) : H^-^ H2, M = H,xH2 

such that 

M{uj) = : xi E Hi}, 

then M{uj) is called a Lipschitz random invariant manifold. If, in addition, for 
every x G H, there exists an x' G M{u) such that for all w G fi, 

\\(j)(t,uj,x) - w,x')||h < ci(x,x', w)e"''2*||x - x'IIh, t > 0, 

where C\ is a positive random variable depending on x and x' , while C2 is a positive 
constant, then M{u) is said to have an exponential tracking property. 

3 Framework 

Consider the following system of stochastic evolutionary equations with two time 
scales 

= -Ax' + -fix', y') + ^w, m H^, (3.1) 
e e Ve 

r = By' + g{x\f), inH^, (3.2) 

where A is a generator of a Co-semigroups on separable Hilbert space Hi, and B 
is a generator of a Cq- groups on separable Hilbert H2. Nonlinearities / and g are 
continuous functions, 

f:HixH2^Hi, g:HixH2^ H^, 

m 

with /(O, 0) = g{0, 0) = 0. The noise process w; = ^ ^j'^jy where {wj}Y=i are two- 

sided Wiener processes taking values in M and hj [1 < j < m) are given elements 
in Hi. Moreover, a is a nonzero constant (noise intensity), and e is a small positive 
parameter representing the ratio of time scales in this fast-slow system. In this 
setting, x' is referred as the "fast" component while y' is the "slow" component. 

Denote by || ■ ||i and || ■ II2 the norms in Hi and H2, respectively. The norm in 
M = Hi X H2 is denoted as || ■ ||. For the linear operators A and B we assume the 
following conditions. 
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(Al) Let A be the generator of a Cq— semigroup e^^ on Hi satisfying 

||e^*x||i < e"^i*||x||i, t>0. 

for all X G Hi, with a constant (i.e., decay rate) 71 > 0. Moreover, B is the 
generator of a Cq— group e^* on H2 satisfying 

||e^*l/||2 < e-^^ii/||2, t < 0. 

for all y E H2, with a constant 72 > 0. 

We also make the following two more assumptions. 

(A2) Lipschitz condition: There exists a positive constant K such that for all 
{xi, Vi) eHiX H2 

- /(x2,?/2)||i < K{\\xi -X2II1 + \\yi -1/2112), 

and 

\\9ixi,yi) - fi'(x2, 2/2)112 < K{\\xi - X2\\i + ||yi - y2\\2)- 

(A3) Assume that the Lipschitz constant K of the nonlinear terms in system 
fl3.1l) - fl3.2l) is smaller than the decay rate 71 of A, that is, 

K < 71. 

Remark 3.1. We note that the system fl3.ip - (l3.2p is an abstract model for vari- 
ous complex systems under random influences, which can be a finite- dimensional, 
stochastic slow-fast system analysed in 13^ \3^. 

Now as in [17], we verify that the stochastic evolutionary system fl3.ip - fl3.2p 
can be transformed into a random evolutionary system which generates a RDS. For 
this purpose, let r/"? be a stationary solution of the linear stochastic evolutionary 
equation 

dr]^(t) = -Ar]^(t)dt + ^dwit). (3.3) 
e y/e 

This means that the random variable 77^ with values in Hi is defined on a {6't}tgig— invariant 
set of full measure such that 

t riH^t(^) 

is a solution version for (13. 3p . Let ^ be the stationary solution of the linear stochas- 
tic evolutionary equation 

d^{t) = A^{t) + adw{t). 

Then by the scale property of Wiener process, r]^{9tuj) has the same distribution 
of ^{6tuj), see the Lemma 3.2 in [32]. For the existence of stationary solutions to 
stochastic evolutionary equations see [5]. 
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Define X'^ = — 77^ {9tuj) and = y'^. Then the original evolutionary system 
fl3.1l) - fl3.2p is converted to the following random evolutionary system 

dX' = -AX'dt + -F(X^ F^ 9'^uj)dt, (3.4) 
dY' = BY'dt + y^ eiuj)dt, (3.5) 

where 

F(x^ r^ ^lo;) = fix^ + T^He^u), r^), 

Let Z^(t,a;,Zo) = {X'{t,uj, Xo,Yo),Y%t,u, Xo,Yo)) be the solution of (M - 
(13.51) with initial data (X''(0), 1^''(0)) = (Xq, Yq) := Zq. Then the solution operator 

of 

$^(t,a;,(Xo,Fo)) = {X^{t,u,Xo,Yo),Y^it,u,Xo,Yo)) 
defines a random dynamical system [TTj. Furthermore 

0'(t,w) := $'(t,w) + (r/^(^ta;),0), t > 0, u eQ 

is the random dynamical system generated by the original system (I3.1l) - (l3.2p . 
We introduce some notations. Let /i be a positive number satisfying 

71 - > K (3.6) 

For any a G M, define Banach spaces 

^a" = \ ^ ■ (~oo,0] I— 7- Hi is continuous and sup ||e~"V(^)lli < 00 

with the norm = sup ||e~"*93(t)||i for i = 1,2. Similarly, we define Banach 

spaces 



= W • [0) 00, ) t-). iJj is continuous and sup ||e ° V^(i)||i < 00 
with the norm ||v2||(^i,+ = sup ||e~"V(^)l|j ^ = Let be the product 



Banach spaces := C^'"^ x C^'^, with the norm 
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4 Exponentially tracking invariant manifolds 



In this section, we prove the existence of a Lipschitz continuous invariant manifolds 
M''{u), with an exponential tracking property, for the random evolutionary system 

Define 

M'{uj) = ^^ZoEM: Z'{-,u,Zo) eCZe} ■ 

This is the set of all initial data through which solutions are bounded by e~^*. We 
shall use Lyapunov- Perron method to prove that M'^{uj) is an invariant manifold 
described by the graph of a Lipschitz function. For this we will need the following 
properties of the random function Z''{-,uj, Zq) (see p!8]). 



Lemma 4.1. Suppose that Z^{-,uj) = cu), cu)) isinC_i,. ThenZ'^{t,u) 

is the solutions of (13.41) - fl3.5p with initial data Zq = {Xq^Yq) if and only if Z''{-,ui) 
satisfies 

\t^e^F{X^{s),Y^{s),9lu)ds \ 




-,Bt 



Yo + /o e^(*-^)G'(X^(s), Y^s), eiu)ds ) 



Theorem 4.1. (Invariant manifolds) 

Assume that (A1)-(A3) hold and that e > zs sufficiently small. Then the random 
dynamical system defined by (I3.4p -f l33|) has a Lipschitz random invariant manifold 
M^(cj) represented as a graph 



where 



is the graph mapping with Lipschitz constant satisfying 

K 



LipH'{uj, ■) < 



(71 - /^) 



1-K 



71-/1 fJ,-e^2 



Proof. We adapt the method of Lyapunov- Perron to fast-slow random dynamical 
systems. To construct an invariant manifold for system (I3.4p - (l3.5p we first consider 
integral equations 




( \r_^e^F{X\s),Y^{s),eiu)ds 
^ e^'Fo + /o e^(*-*)G(X^(s), Y\s), dlu)ds 



, t < 0. (4.1) 



A solution of (gl]) is denoted by Z'{t,uj,Zo) = {X'{t,u,Yo),Y'{t,u,Yo)) . Intro- 
duce the operators JT"/ : C"m C^'^ and ^2 ■ CZe ^ by means of 

^ J — oo 

Jo 

for t < and define the mapping J''^ by 



It can be verified that J'^ maps C_m into itself. To this end, taking z = {x,y) G 
C~M, we have that 



e t<o 



oo 



< -supI / e(^+^)(*-)rf4lkllc- 
e t<o W_oo j -# 

kllc_^, (4.2) 
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and 



-T t<0 { Jt ) 

< i^supj f e(-^2+^)(*"")dsl||z||^- 

t<o y Jt } ~ 



\c- 



A\c-.r (4-3) 



/i-e72 

Hence, by definition of J'^ we obtain 

\\J\z)\\c-^ <K(ir,7i,72,/i,e)||2;||^-^ 

with 

K eK 

/t(i^,7i,72,/i,e) = \ . 

7i - /i /X - e72 

Thus, we conclude that J'^ maps C~£, into itself. 



Next we show that the mapping J*" is contractive. To this end, taking z 
(x, y),z = (x, y) G C"^ , we have that 



+ ||l/(s)-y(s)||2)rfs 

< -sup I r e("?"+^)(*-)d4lk- 



^ >-^llc_,, (4.4) 



7i -/i 



and 



11^2(2) -^2 (^)llc;2.- < Ksuple^' [ e-^'^'-'^e-^'ds\\\z - z 

~T t<0 I Jt J 

t<o I it J 



IC" 



^ ^ ^\\c- 



eK „ „ , , 

z-z\\c-^- (4-5) 



/i - £72 
Hence, by (lOj) and f l431) 

\\T{z) - Tmc:^ < <K,^,,^,,fi,e)\\z - z\\c-^, 

where 

K(ir,7i,72,/i, e) = ^ > 



7i - /i /i - 672 7i - 

as e — )■ 0. Taking into account of (13. 6p there is a sufficiently small constant eo > 
such that 

k{K, 7i, 72, /i, e) < 1, /or e G (0, eo]. 
Therefore, the mapping J''^ is strictly contractive in C"^, , and, consequently, the in- 

tegral equation (14. ip has a unique solution Yo) = {X'{t, u, Yq), Y'{t, u, Fq)) 

in C~u. Furthermore one has the estimate 



\Z\,u;,Y,)-Z^{.,u;,Y,)\\c-^ <- 7- 7^1-1^2112 (4.6) 

1 - k(K, 71,72, /i,e) 



for all CO e n,Yi,Y2 E H2 
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Define 

H'{u,Yo) = - f e-^'/'F{X'{s,u,Yo),Y'{s,u,Yo),9lu)ds, (4.7) 

we then get from (14. 6p 

\\H^iu,Y,)-H^iu,Y2)h< ( ^ . T, TjT^ rA\yi-y2h 

for all ^1,^2 & H2,UJ E Q. It then follows from Lemma [4.11 that 

M'{u) = {{H'{u,Yo),Yo) : Fq G H2} . 
In order to see that M^{u) is a random set we need to show that for any z = 

u ^ mf^||(x,y) - (4-8) 

is measurable, see Castaing and Valadier [7|, Theorem III. 9. Let He be a countable 
dense set of the separable space H. Then the right hand side of (14. 8 p is equal to 

inf \\ix,y)-mu,Vz'),Vz')\\ (4.9) 

which follows immediately by the continuity of iJ^(a;, ■). The measurability of any 
expression under the infimum of (14. 8 p follows since u — )■ H'^{u,Vz') is measurable 
for any z' G H. 

It remains to show that M'^{u) is invariant, i.e., for each Zq = {Xq, Yq) G M'^{u), 
Z%s,uj,Zo) E M^eicu) for all s > 0. We first note that for each fixed s > 0, 
Z'^{t + s, u, Zq) is a solution of 

dX' = -AX'dt + -F(X^ y^ 9U9':u))dt, 
e e 

dY' = BY'dt + G'(X^ ei{9lu))dt, 

with initial datum Z(0) = {X{0),Y{0)) = Z'{s,u,Zo). Thus, Z'{t + s,uj,Zo) = 
Z%t,eiLO,Z\s,cu,Zo)). Since Z%-,co,Zo) e CZ,, Z^t^Olu, Z%s,u, Zq)) G CZ,. 

e e 

Therefore, Z''{s,uj, Zq) G M''{61uj). This completes the proof. □ 



Remark 4.1. We remark that the key point in the proof of Theorem 4-i is that 
K(K,7i,72,/i,e) = \ < 1. 

7l - yU - 672 

In the particular case where e = 1, one has k = — — I — — < 1, which is the usual 
spectral gap condition. We note also that the proof is valid for sufficiently small 
e > only in the case t^^zji < 1- This explains the assumption (A3). It is unclear 
to us about how to relax this condition. 
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In what follows we prove the exponential tracking property which means the 
manifold M'^(a;) attracts exponentially all the orbits of $^ on condition that the 
scaling parameter is sufficiently small. 



Theorem 4.2. (Exponential tracking property) 

Assume that the assumptions (A1)-(A3) hold. Then for sufficiently small e > 0, 
the Lipschitz invariant manifold for (I3.4p - fl3.5l) obtained in Theorem 4-i ho,s the 
exponential tracking property in the following sense: There exist constants Ci > 
and 6*2 > such that for any Zq = (Xq, Yq) G H there is a Zq = {Xq, Yq) G M''{u) 
such that 

\\^'{t,uj, Zo) - ^'{t,uj, Zo)\\ < Cie-^^iZo - Zo||, t > 0, 

where || ■ || denotes the norm in space M. = Hi x H2 defined by 

\\4 = + Wvh, z = {x,y). 

Proof Assume that Z'{t) = {X'{t),Y'{t)) and Z'{t) = {X'{t),Y'{t)) are two 
solutions for (^-(^, then Z'{t) = Z'{t) - Z'{t) := {U'{t),V'{t)) satisfies the 
equations 

dU' = ^AU'dt + ^F{U% V\ 9lu)dt, (4.10) 
dV = BV'dt + G(f/^ 1/^ eiuj)dt, (4.11) 

where 

F{u% v% ei) = F{u' + x\ v + y^ eiu) - f{x\ y\ eiu), 

and 

G{u', v\ ei) = G{u' + x% v + y^ eico) - G(x^ y^ eico). 

First we claim that Z'{t) = {U'{t), V'{t)) is a solution of flCTl - fHTTl) in C^^ if 




e^*A[/^(0) + i e^(*-*)F(t/^(s), V'{s), 9lu)ds 
Jl e^(*-^)G(f/^(s), V'{s), eiu)ds 



(4.12) 



This can be verified by using the variation of constants formula. Next we are 
going to prove that Km has solutions {U',V') in C^^^ with (f/^(0), 1/^(0)) = 

(?7(0),y(0)) and such that (Xq, Fq) = (f/(0),y(0)) + (Xo,Fo) e M'{u). Recall 
that 

- I 

(Xo,?o) e M\uj) ^Xo = - e''^-'^F{X'{s,Yo),Y'is,Yo),eiu)ds. 

^ J —00 
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It follows that 

(Xo,Fo) = (f/(0), V^(O)) + (Xo,Fo) e M^uj) 

if and only if 

1 

U{Q) = -Xo + - e''^~''^F{X'{s,ViO) + Yo),Y^is,ViO) + Yo),eiu)ds 
e J-oo 

= -Xo + H%u,ViO) + Yo). (4.13) 
For every Z = {U, V) G define for t > 

1 /•* 

Xi^(Z(-))[t] := e^*/^f/(0) + - / e^(*-^)/^F(t/(s), 

^ Jo 

and 

X^(Z(-))[t] := f e^^'-^^'^G{U{s),V{s),eiuj)ds, 

J oo 

where t/(0) is given by fl4.13l) . Define the operator X*^ by 



T{Z{.)) 



^1{Z{-)) 



It is easy to see that X^ is well-defined from to itself. To this end, assume 
that Z^Z ^ C'^ii; we obtain from f l4.13p the estimate 

e 

||e^*/^(f/(0) - f7(0))||i < e-^^^l'UvE%V{^)-Y{^)\\2 



j e^(-^) {g{Z{s), 9lu) - G{Z{s), Olu)) ds 

J oo 



' oo 

r"00 



POO 

< e-^^^''LipH' ■ K I e'^'WZis) - Z{s)\\ds, 
Jo 



and so 



\Il{Z - Z)V.+,i < UpH' ■ K\\Z - Z\\c+ sup(e-(-^+^)* /"e(^^-^)^rfsj 

~T -T t>0 [Jo J 



+ ^||Z-Z||^+^sup<|e^* / e-^'^'-''>/^e-T^ds 



< (—^ + ]\\Z-Z\\c+ . 4.14 

For the operator Xg we have 

||X^(Z-Z)||^+,2 < A'llZ- sup|e-(-^+T2)* /" e(-^+^2)^cisj 
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< 



eK 
/^-72 



Recalling that 



LipH' < 



K 



(4.15) 



(71 - /^) 



and taking fl4.14p and fl4.15p into account, we obtain 



\XHZ - Z) 



< p(ir,7i,72,/i,e)||Z-Z| 



with 



p{K, 7i,72,/i,e) 



K eK 

+ 



7i - /i /i - 672 



(71 - /^) (7 - 72) 
K 



l-K 



+ 



71— /X /I— €72 



— 



7i -/i 



as e — )■ 0. By fl3.6p there is a sufficiently small constant e^, > such that 
p(_ft', 7i, 72, /i, e) < 1 for all < e < eg. Therefore, the operator is strictly 
contractive and has a unique fixed point Z G C^m which is the unique solution 

for fimi) and satisfies (Xo,>o) = (t/(0), 1/(0)) + (Xo, >o) e M^(u;). Moreover, we 
have 



12^1 



< 



1 



1 



eK 



which means 

||$^(t,u;,Zo)-$^(t,a;,Zo)|| < 



71—/^ M-e72 



e £ 



1^(0)11 



_|_ eif 



l^o-Zoll, t> 0. 



^7l-/i M-e72_ 

Therefore, the exponential tracking property of M^(ci;) is obtained. 



□ 



Remark 4.2. 5y the relationship between solutions of system fl3.ip - fl3.2p and 
fl3.4p - fl3.5p . i/ie original fast-slow stochastic system also has a Lipschitz random 



invariant manifold under the conditions of Theorem\4. 1\ which is represented as 



M' ^ 



UJ 



.0) 



H2] 
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with 

h'{u,Yo) = H^iu,Yo)+r]7(co). 
Hence, if system (13. 4p - fl3.5p has an exponential tracking manifold so has system 



Remark 4.3. Theorem implies that any orbit of the fast-slow system tends 
exponentially to an orbit on the manifold M^(a;) which is governed by an evolu- 
tionary equation with usual time scale. To be more specific, i.e., for any solution 
= {X',Y') for there is an orbit Z%t,u) = {X'{t,u),Y%t,u)) on 

the manifold M*" which satisfies the evolutionary equation 



Y' = BY' + G [h^OIu, Y'), Y\ eiu 



such that 



— tLf 

e f 



Z'{t,u)- Z'{t,uj)\\ < ^||Zo-Zo||, t>0 



K _|_ eK 
71 -A* M-e72 

where Zq = (X^(0), F^(0)) and Zq = (X(0),y(0)). 

Applying the ideas from the Remark 14.21 we have a reduction system which 
describes the long-time behavior for system fl3.ip - fl3.2l) . 

Theorem 4.3. (Reduction system) 

Assume that e > is sufficiently small and the assumption (Al )-(A3) hold. Then 
for any solution z^{t) = {x^ (t) , (t)) with initial data z'^{0) = (xo,yo) to system 
(O-dS^), there exists a solution z'{t) = {x'{t),y'{t)) with z{0) = (F(0),y'(0)) = 
(xo,yo) to the reduced system 

y' = W + 9 {x, 
x = h%eiu,y^), 

such that for any t > and almost sure u E Q, 

\\z\t,u:) - z\t,uj)\\ < ^— —^\\zo-Zo\\ 

1 _ (_f^ + ) 

— fit ~ I I 

< CK,'yi,fj,e~\\Zo — ZqW 

with — < and Ck,'^^,!! being a constant depending on K,'ji and /i. 
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5 Slow manifolds 

Now we consider an asymptotic approximation for the invariant manifold M'^{u), 
as e — 0. 

The scahng t — )■ et in system f l3.4p -( l375ll yields 

dX' = AX'dt + Y% eituj)dt, (5.1) 

dY' = eBY'dt + eG{X% Y% eitUj)dt, (5.2) 

where 

F(X^F^e» = f{X^ + r]7(e,tu),Y^), 
G{X% Y\ ei,uj) = g{X^ + /^'(^eiu;), Y^). 

We now replace rj^{9^iu) by ^{OtOj) that has the same distribution, then we have a 
random evolutionary system, whose solution's distribution coincides with that of 
the system ( 15.ip -( !5^ . in the form of 

dX' = AX'dt + f{X' + ^{Otuj), Y')dt, (5.3) 
dY' = eBY'dt + eg{X' + ^(^tw), Y')dt. (5.4) 

By proceeding as in the proof of Theorem 14.11 it can be shown that the system 
( ]5.3p -( !5^ has a random invariant manifold represented as 



with 



where 



H'{u,Yo)= f e^'f{X^{s,uj,Yo)+i{etu:),Y\s,u,Yo))ds, 



X\t,u,Y^)= f e^'^'-'^f{X%s,u,Yo)+i{e,u),Y^{s,u,Yo))ds,t<Q, 
J —00 

Y\t, u, Yo) = e^'^Fo + e / e^(*-^)^(?(l^(s, u, Yo) + m^), Y^s, u, Yo))ds, t < 0, 

is the unique solution in CZ^ for the above integral equations. With a change of 
variables s/e — t in (14. 7p , we have 



H\u,Yo) = I e-^J{X'{se,co,Yo) + vHdetco),y'ise,u,Yo))ds 
J —00 

'f{X'{se, 00, Yo) + vHOetOo), Y'{se, u, Yo))ds, 



-00 
—00 
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where ^ denotes equivalence (coincidence) in distribution. Therefore, the invariant 
manifold M^{u) is a version in distribution for M'^{u). 

Next, we show that there exists a random invariant manifold M^{u!), which is 
called a random slow manifold for system fl5.3p - fl5.4l) . will be the asymptotic limit 
of the manifold M^(cj) as e — )■ 0. To this end, we consider the following system 

dX = AXdt + f{X + ^{etUj),Y)dt, (5.5) 
dY = 0. (5.6) 

By the same discussion in Theorem 14.11 system fl5.5l) - fl5.6p has a random invariant 
manifold with representation 

M\u) = { {H\u, Yo), Yo) ■.YoeH2}, (5.7) 

where 

H\uj,Yo) = / e-^'f{X{s,uj,Yo)+adsOo),Yo)ds, 

J — oo 

and X{t,u,Yo) is the unique solution in Ci'~ for integral equation 

X{t,u,Yo)= I e^^'-'^f{X{s,uj,Yo)+adsUj),Yo)ds, t<0. 

The main result of this section is the following theorem. 
Theorem 5.1. (Slow manifolds) 

Let the assumptions (A1)-(A3) hold and also assume that there exists a positive 
number Cg such that sup ||(7(x, ?/) ||h2 = Cg. Then the invariant manifold 

M'^{uj) for the system fl5.1|] - fl5.2p can be approximated by a slow manifold M^{u) 
defined in {\5.7\i . in the sense that their respective graph mappings H"^ and H satisfy 

\\H^{u,Yo)-H{u,Yo)h = 0{e), 

or 

H'{u,Yo) = H{u,Yo) + 0{e), 
for all Yq G V{B), a.s. u ^ Q and as e — t- 0. 

Proof. In this proof, the letter C with or without subscripts denotes positive con- 
stants whose value may change in different occasions. We will write the depen- 
dence of constant on parameters explicitly if it is essential. As is known [28], if 
YoeV{B),t<0, 

\l„Bte^^ II _ II / R 



Yo-Yoh = II / e^Wodrjb 

Jet 
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Jet 



||S>o||2-(e-^^^*-l). (5i 

72 



Then we have, for alH < 0, 
\\Y'{t,oo,Yo)-Yo\\2 < ||e^^*Fo-Fo||2 



l) + eCgj e-'''^^'-'^ds 



1 

^"0112- 
72 



< \\BY.\\o-ie-^^'' 



= C(e^^2^*-1). (5.9) 
Then, by using again (15. Sp . we have 

\\X^{t,oo,Yo)~X{t,uj,Yo)\\i < K f e-^^(*-^)||X^(s,a;,Fo)-X(s,a;,Fo)||irfs 

J — oo 

+KC I e-"<'^'-'\e-^^'' - l)ds 
= K f e-'^'-'-''^\\X^{s,uj,Yo)-X{s,uj,Yo)\\ids 

J — oo 



+C{ i e-'^'' - — 

7i - £72 7i ■ 



which imphes 



\\X^{-,U,Yo)-X{;U,Yo)\U- < K\\X^{;UJ,Yo)-X{;U,Yo)\U 



sup / e 

t<0 



oo 



+Csup{e'^*( i e-'^'' - — )| 

t<o 7i - £72 7i 

— ^||1^(-, a;, Yo) - X{; u, Yo) 11^1,- + C sup ^(t, e), 
7i - At t<o 

(5.10) 



where 



y{t,e) = e^*( - — ), t G (-oo,0], 

7i - £72 7i 



with 

^(0,e) = ^ 1. (5.11) 

7i - £72 7i 
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Furthermore, for sufficiently small e > 0, 

d-^jt, e) ^ ^^f/ ^-e72 ^-e^,t /i x 
dt Si -£72 7i 

7i-e72 7i 
^^^ -€72(71 -^) 
71(71 - £72) 

< 0, t e (-cx),0). (5.12) 
Now, according to flS.lOp . f lS.lip and fl5.12p . we have 

\\X^{;u,Yo)-X{;u,Yo)\y,- < -^\\X^ {■ , CO ,Yo) - X {■ , u ,Yo)\y.- 

7i - e72 7i 

By iH), 

||l^(-,a;,Fo)-X(-,a;,Fo)|lci.- < C( ^ -). (5.13) 

7i - £72 7i 



Hence, thanks to (15.91) and (15.131) . we deduce 



r-O 

\H^{u;,Yo)-H{co,Yo)h < Ki e^''\\X%s,uj,Y,) - X{s,uj,Y,)\\,ds 



00 



+K I e<''\\Y\s,u,YQ)-Y42ds 
J —00 

< C( ^ — ) r e^-^'-^^'ds 

7i - £72 7i i-00 

+C /" e^i^(e-^2" - l)ds 
J —00 

= ^^(^ -) = 0{e). 

7i - £72 7i 

This completes the proof. □ 

Remark 5.1. Consider the case that H2 is a finite dimensional space, the operator 
B is a constant matrix and (I5.3p -f l5^ is a coupled system of an evolutionary 
equation and ordinary differential equations. This system arises from biology, such 
as Hodgkin- Huxley systems {see Example {KM- Then the above theorem implies 
that for any bounded set E C H2, 

sup \\H\u,Yq) - H{uj,Yq)\\i = C(e), a.s. ueQase^O. 

YoGE 
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6 Illustrative examples 



Let us look at several examples to illustrate the results in the previous two sections. 

Example 6.1. Let D G M.^ be a bounded domain with smooth boundary dD. Con- 
sider a coupled system of stochastic parabolic-hyperbolic partial differential equa- 
tions {see, e.g., [6] and [16] ) 



111 
ut = -(Am - au) + -f{u, v, Vt) + —=w{t), (6.1) 
e e y/e 

u = on dD, (6.2) 

vtt = - (3v + g{u,v,vt), (6.3) 

v = on dD, (6.4) 

where A denotes the Laplace operator and the parameters a, /3 are positive. The 
interaction functions 

/ : ^ M and g : ^ R 

are assumed to be Lipschitz continuous with a Lipschitz constant K > 0. Thus 
the assumption (A2) holds. Such a system may describe a thermoelastic wave 
propagation in a random medium fT^ . The wave profile v in an interacting random 
thermoelastic medium is described by a hyperbolic partial differential equation. If 
the wave is temperature dependent and the heat conductivity has faster evolution, 
then the hyperbolic equation is coupled to a stochastic parabolic (heat) equation with 
different characteristic timescales. 

We introduce the usual Hilbert space L'^{D) as well as the Sobolev spaces H^{D) 
and Hl{D). Take Hi = L'^{D). Let A = A-aha with domain V{A) = H'^DH^. By 
the semigroup theory the operator A generates a contraction semigroup je"^* : t > 
0} inHi ([21]) which satisfies \\e^^\\Hi < e~^^\t> with-fi = a. Let B = A-/3Iid 
with domain D{B) = H"^ fl H^. Define 



B 

and H2 = Hq{D) x L'^{D) with the energy norm 

\\4h2 = {Mm + Wv'Wh} 



where \\ ■ ||^i and \\ ■ \\l2 denote the norm in Hq and L^, respectively. Let T>{E) = 
'D{B) X H^. It is known that B generates a unitary group ([33]) in H2 which satisfies 
\\^'^'^\\h2 ^ e"'^2* fort G M with 72 = 0. Then the system fl6.ip - fl6.4p can be rewritten 
as 

Ut = -Au + f{u, z) + -^wt, 
e ye 
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Zt = Bz + G{u, z), 

with 

G{u,z) = (0,g{u,z)), 
which is in the standard form of (13. ip - fl3.2p . Thus under the condition 

K < 71 

and the scaling parameter e small enough the random dynamical system generated 
by fl6.ip - fl6.4p has an invariant manifold Ai'^{u) = {{h'^{uj,YQ),YQ) Yq G H2}, pos- 
sessing the exponential tracking property by Theorem \4.S\ Moreover, by Theorem 
\4-3\ the reduction system for long-time behavior to system f l6.ip - fl6.4p is 

t = By^ + Gix,y^), 
x = h%eico,y^). 

Note that a similar result for this example has also been obtained in [B]. 

Example 6.2. Let D G be a bounded domain with smooth boundary dD. 
Consider a system of coupled parabolic partial differential equations and ordinary 
differential equations 




u = on dD, (6.6) 
vt = g{u,v), (6.7) 

where f : M}~^"^ R, g : M}^^'"- 1— )■ are Lipschitz maps with a Lipschitz constant 
K > 0: 

\fixi,yi) - /(x2, 1/2)1 < K{\xi - X2I + \yi - 2/2k'"), 
\g{xi,yi) - 5' (2:2, 1/2) < Ki\xi - X2I + \yi - y2\m--), 

for all {x,y) G M x M™. Thus the assumption (A2) holds. This system may model 
certain biological processes. For instance, the famous FitzHugh-N agumo system 
\2^ , as a simplified version of the Hodgkin-Huxley model fT^ - which describes 
mechanisms of a neural excitability and excitation for macro-receptors, belongs to 
this class. 

As in Example \6.1\ the differential operator A = A with domain V{A) = H'^dHQ 
generates a Co-semigroup {e^* : ^ > 0} on Hi = L'^{D) which satisfies ||e"^*||_H'i < 
e-'^i* with 71 = Mspec{A} > 0. Let B = in H2 = [L'^{D)]"'. It is clear that 
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e^* = lid for a// t G M and ||e^*||i^2 ^ e"'*'^* with 72 = 0. Therefore, the system 
(16.51) - fl6.7p has a random invariant manifold Ai'^{uj) = {(/i'^(a;, lo); ^0) ^0 ^ H2} 
with an exponential tracking property if K < '-fi and e > is small enough. We 
also have the reduction system 

x = h^{9lco,y^), 
for the long time behavior of the original system (16. 51) - (16. 71) . 



Example 6.3. Consider the following system of two coupled wave equations (i.e., 
hyperbolic partial differential equations) on a bounded spatial interval I = [0, tt] : 

111 

utt = -{Au - vut) + -/(u, vt) + -^w{t), (6.8) 
e e a/ e 

u = on dl, (6.9) 

vtt = Av - (3v + g{u,v,vt), (6.10) 

v = on dl, (6.11) 

where A denotes the Laplace operator and the parameters (3, v are positive. The 
interaction functions 

/ : ^ M and ^ : M 

are Lipschitz continuous with a Lipschitz constant K > 0. Thus the assumption 
(A2) holds. This system models, for example, vibrating strings connected in par- 
allel with zero boundary conditions IMf and multi- component wave fields such as 
electromagnetic waves in plasmas, elastic waves in solids, light waves in anisotropic 
and inhomogeneous media f2B^ . 

Rewrite the equations (I6.8p -( l6l9|) as 

^ = lA^U + -F{U,V) + ^W{t), 
at e e y/e 



where 



and 



el,d \ ... f \ / 



U = {u, u'),V= {v, v') e iJo (0, tt) X L\0, vr). 
The linear operator A"^ has the eigenvalues 
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with the corresponding eigenvectors 




sin kx 
sin kx 



) 



It is clear that the operator generates a Cq— semigroup e-^'* on Hilhert space 
Hi := Hq{Q,tt) X L^(0,7r) equipped with energy norm introduced in Example \6.1\ 
and it satisfies 

We^^'Wn, < t > 0. 

In the same way as in Example \ 6.1\ the linear part of the equation f l6.10p -f l6.lip 
generates a unitary Cq— semigroup on Hilhert space H2 = Hq{0, vr) x L^(0, vr). Thus 
under the condition that K < u, the system fl6.8p - fl6.lll) has an exponentially 
tracking random invariant manifold A4.'^{u) = {{h'^{u, Yq), Yq) Yq E H2} when e > 
is sufficiently small. In particular, by Theorem \4.3\ the system f l6.8p -f ETTT]) has a 
reduction equation 



where B and G are defined as in Example \6.1[ 

7 Remarks on the case of local Lipschitz nonlin- 
earity 

We have limited ourselves to the case where the nonlinearities are globally Lips- 
chitz continuous. We remark that when the nonlinearities in those three examples 
in Section 6 are only locally Lipschitz (say, near the origin (0,0)), the above dis- 
cussions remain valid locally. To this end we state the definition of a local random 
invariant manifold [H |8] . 

Definition 7.1. We say that the random dynamical system (j){t,uj) has a local ran- 
dom invariant manifold (LRIM) with radius R, if there is a random set M.^{uj), 
which is defined by the graph of a random continuous function ip{oj, ■) : -B_r(0) f]H2 ^ 
Hi, such that for all bounded sets B in Bfj{0) C H2 we have 




(P{t,u)[M\u)[]B] c M^'iOtu) 



for all t G (0, Tq{uj)) with 



Tq{uj) = ro(a;,5) =inf{t > : 0(t,a;)[>l V) f|5] /Sr(0)}. 
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Let X '■ Hi X H2 ^ he a bounded smooth function such that 

X(^l,^2) = 



if + ||t;2||2 < 1, 
if ||fi||i + balh > 2. 



For any positive parameter R, we define Xr{vi,V2) = xi'^j^) for (^'1)^2) G 
Hi X H2. Let f^^\x,y) := XRix,y)f{x,y),g'^^\x,y) := XB{x,y)g{x,y). For every 
R> 0, there must exist a positive Kji such that 

- f''\x2,y2)\\i < KniWxi - x^Hi + \\yi - 2/2II2), 

and 

2/1) (0:2,2/2)112 < Kr{\\xi - X2II1 + Ibi - y2||2). 
Then the cut-off system of fl3.4p - fl3.5p is as follows: 

dX' = ^AX'dt+^F^^\X%Y',9'tu)dt, (7.1) 

dY' = BY'dt + G^^) (X^ F^ ^,^a;)cit, (7.2) 

where 

F(^)(X^r^^,^a;) = /(^)(X^ + r/^(M,>^'), 
G'^^\X',Y',etuj) = g'^^\X' + 7]^etUj),Y'). 

The system f l7.ip -( rr!2l) has a unique solution and thus the solution mapping gen- 
erates a continuous random dynamical system If Kji < 71, then the cut-off 
system f l7.ip - fl7.2p admits a globally invariant manifold M% possessing exponen- 
tially tracking property. Now as and $^ agree on i?/j(0), we conclude that 
= n Hr{0) defines a local invariant manifold of the original system (13. 4p - 
(O. 
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